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I  INTRODUCTION  AND  SUMMARY 


Recent  studies  using  the  DNA  Wideband  satellite  have  indicated  that 

the  transionospheric  radio  channel  is  best  represented  by  a  two -component 

multiplicative  model  that  separates  received  signal  scintillations  on  the 

basis  of  their  time  structure  (Fremouw  et  al.,  1976,  1978).  This  model 

appears  to  be  useful  from  VHF  to  L-Band,  and  is  sufficiently  general  that 

many  simple  channel  representations  (e.g.,  Rayleigh,  Rician,  log-normal) 

can  be  obtained  as  special  cases.  In  this  report  we  develop  a  general 

expression,  applicable  to  a  large  class  of  binary  diversity  receivers, 

for  the  probability  of  bit  error  (P  )  using  this  channel.  We  then  special- 

£ 

ize  this  expression  for  noncoherent  FSK  and  dif ferentially-coherent  phase- 
shift-keying  (DPSK)  communication  systems,  and  examine  the  behavior  of  P£ 
as  the  scintillation  parameters  are  varied.  The  expressions  for  Pg  were 
evaluated  numerically  on  a  CDC  6400  computer;  closed-form  expressions  for 
Pg  could  be  obtained  only  in  the  simplest  cases. 

The  results  all  tend  to  show  that  amplitude  scintillation  (fading) 

is  by  far  the  dominant  factor  in  the  determination  of  the  error  rate. 

Phase  scintillation  has  some  effect  on  DPSK  if  it  is  rapid,  but  it  has 

virtually  no  effect  on  either  FSK  or  DPSK  if  it  is  slow.  In  line  with 

this,  we  found  that,  for  reasonable  parameter  values,  the  focus  component 

has  relatively  little  effect  on  P  ,  because  it  introduces  primarily  slow- 

e 

phase  scintillation.  The  addition  of  a  coherent  component  to  a  scintil¬ 
lating  signal  reduces  the  error  rate,  of  course,  but  significant  reduc¬ 
tions  are  not  obtained  until  the  coherent  component  contains  more  than 
one-half  the  total  signal  energy. 

As  with  Rayleigh  fading,  rapid  scintillation  always  increases  the 
error  rate.  Rician  and  non-Rician  models  show  essentially  the  same  de¬ 
pendence  on  the  fading  rate  (with  one  exception)  as  Rayleigh  fading.  The 


A  list  of  references  will  be  found  at  the  end  of  the  report. 


exception  occurs  when  most  of  the  signal  energy  is  in  the  quadrature 
term  of  the  scatter  component.  Under  these  conditions,  DPSK  becomes 
much  more  sensitive  to  rapid  scintillation  than  FSK. 

Channel  Model 

The  channel  model  proposed  by  Fremouw  et  al.  (1976a, b)  represents 
the  complex  received  signal  in  the  form 

r(t)  =  rs(t)  rf(t)s(t)  +  n(t)  ,  (1) 

where  s(t)  is  the  transmitted  signal,  n(t)  is  additive  noise,  and  rg(t) 

and  r^(t)  are  two  statistically  independent,  complex  stochastic  processes 

having  Gaussian  and  log-normal  statistics  respectively.  In  the  transiono- 

spheric  channel,  r  (t)  represents  the  diffractive  scattering  from  small- 
s 

scale  irregularities  and  is  referred  to  as  the  "scatter"  component.  The 
log-normal  term,  r^(t)  represents  the  refractive  focusing  from  large 
scale  irregularities,  and  is  referred  to  as  the  "focus"  component.  The 
correlation  time  for  the  scatter  component  is  much  smaller  than  the  cor¬ 
relation  time  for  r^(t);  hence,  the  scatter  term  is  associated  with  "fast 
fading,"  or  rapid  scintillation,  and  the  focus  term  is  associated  with 
slow  variations  in  the  average  amplitude  and  phase.  A  functional  model 
for  the  channel  consists  of  two  cascaded  multiplicative  filters  as  shown 
in  Figure  1. 

The  scatter  component  is  given  by 

rg(t)  =  T]  +  hs(t;f)  ,  (2) 

where  7]  is  a  coherent  or  specular  term,  and 

hg(t;f)  =  xs(t;f)  +  jys(t;f)  .  (3) 

The  quadrature  terms  xg(t;f)  and  yg(t;f)  are  zero-mean,  stationary  joint 
Gaussian  processes  with  variances 
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FIGURE  1  FUNCTIONAL  DIAGRAM  OF 
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THE  TWO-COMPONENT  MULTIPLICATIVE 


=  a 


=  CT 


(4) 


and  correlation  coefficient,  p  ,  given  by 


p  =  x  (t;f)  y  C t ; f )/o  a 
xy  s  s  x  y 


(5) 


The  overbar  indicates  the  average  (ensemble)  of  the  quantity  beneath.  It 
is  assumed  that  the  functions  involved  are  ergodic.  The  focus  component 
is  given  by 


r£(t)  -  e*(t'£) 


(6) 


where 


ty(t;f )  =  X(t;f)  +  j$(t;f)  +  e 
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In  Eq.  (6),  x(t;f)  and  $(t;f)  are  zero-mean  stationary  joint  Gaussian 

2  2 

processes  with  variances  ct.  and  cr  and  correlation  coefficient  p 

X  $  X$ 

We  make  the  simplifying  assumption  that  hg(t;f)  and  ij<(t;f)  are 

constant  over  the  frequency  band  of  interest,  and  shall  henceforth  drop 

the  f  variable  from  our  notation.  Wi  also  assume  that  the  channel  is 

lossless. 

This  imposes  the  constraints 


T]2  +  |hs(t)|2  =  l]2  +  a2  =  1 


(7) 


2  2  2 

where  a,  =  ct  +  a  ,  and 
h  x  y 


,t(t) 


=  1 


(8) 


It  is  straightforward  to  show  that  Eq.  (8)  implies  that 


e  =  -  X  = 


(9) 


Seven  parameters  are  required  to  characterize  this  channel: 

2  2 

•  ct  ,  a  and  p  are  the  variances  and  correlation  coefficient 

x  y  xy 

of  the  scatter  component, 

2  2 

•  ct  ,  ct  and  p  are  the  variances  and  correlation  coefficient 

X  •  X« 

of  the  focus  component,  and 

•  T  is  the  correlation  time  of  the  scatter  component. 

2 

The  coherent  intensity  T|  can  be  determined  from  (7).  The  properties  of 

2  2 

the  scatter  and  focus  components  are  such  that  ct  «  a  .  and  the  corre- 

X  $ 

lation  time  of  the  focus  component  is  much  larger  than  that  of  the  scatter 
component.  Thus,  for  the  bit  error  rate  computations  only  t  need  be 
specified. 


A  typical  example  of  a  signal  undergoing  weak  scintillation  is  shown 
in  Figure  2.  The  scatter  and  focus  components  are  shown  separately  in 
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Figures  3  and  4  (see  Fremouw  et  al.,  1976a,  for  details  of  the  separation 

procedure).  From  the  scatter  diagram  shown  in  the  top  part  of  Figure  3 
2  2 

we  see  that  a,  <  T|  .  Moreover,  because  of  the  rotated  elliptical  form 
h  2  2 

of  the  scatter  diagram  a  <  a  and  the  correlation  coefficient,  p  ,  is 

x  y  xy 

negative  (Rino  and  Fremouw,  1973;  Rino  et  al. ,  1976). 


From  the  scatter  diagram  in  the  top  part  of  Figure  4,  we  see  simi- 
2  2 

larly  that  «  a ^ .  A  careful  inspection  will  show  that  the  amplitude 
and  phase  are  anticorrelated.  From  a  ray-optics  viewpoint,  this  is  ex¬ 
pected  since  phase  depletions  cause  focusing  and  vice  versa.  Indeed, 
the  same  argument  can  be  applied  to  reconcile  the  fact  that  p  <  0; 
since,  for  weak  scatter  (cr^  «  1),  the  phase-quadrature  component  approxi¬ 
mates  the  signal  phase  perturbation,  and  the  in-phase  component  approxi¬ 
mates  the  fractional  amplitude  perturbation. 


An  example  of  severe  scintillation  is  shown  in  Figure  5.  The  scatter 
diagram  for  the  composite  signal  [Figure  5(a)]  appears  to  be  Rayleigh.  We 
see  from  Figures  5(b)  and  5(c),  however,  that  the  scatter  component  re¬ 
tains  a  measurable  coherent  part,  and  the  focus  component  is  dominated 
by  large  slow  phase  variations. 


Thus,  while  the  overall  fading  structure  of  the  channel  for  severe 
scintillation  is  reasonably  well  approximated  by  a  Rayleigh  distribution, 
the  time  structure  of  the  signal  amplitude  and  phase  are  very  different. 
The  analysis  presented  in  the  remainder  of  this  section  lays  the  founda¬ 
tions  for  calculating  the  impact  of  this  model  on  representative  communi¬ 
cation  systems. 
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PHASE  AMPLITUDE  Im  PART  Re  PART 


II  ERROR -RATE  DETERMINATION 


Derivation 

In  this  section,  we  shall  obtain  an  expression  for  the  probability 
of  error  over  the  two-component  channel  for  the  canonic  receiver  shown 
in  Figure  6.  This  canonic  receiver  is  a  general  structure  that  can  be 
made  mathematically  equivalent  to  several  types  of  binary  diversity  re¬ 
ceivers  (such  as  FSK,  DPSK,  fixed  reference  PSK)  by  the  proper  choice  of 
the  two  input  filters  and  the  quadratic  combiner  (see  Appendix  A).  The 
method  is  a  generalization  of  the  procedure  developed  by  Bello  and  Nelin 
(1962a, b).  We  assume  that  the  outputs  of  the  diversity  branches,  d^, 
are  statistically  independent  and  identically  distributed.  Complex 
signal  notation  is  used  throughout. 


INFORMATION 

FILTER 


REFERENCE 

FILTER 


FIGURE  6  CANONIC  BINARY  RECEIVER  STRUCTURE 


The  outputs  u^  and  v  ,  of  the  information  and  reference  filters  in 
Figure  6  can  be  expressed  as 


Uk  =/r k(C)  8l(t) 


dt 


(10) 


and 


b 


vk =  ;rk(t)  go(t)  dt 


(in 


where  the  Limits  of  integration  depend  upon  and  gQ  (see  appendix). 
If  we  make  the  definitions 


\  =  Ukl  +  jUk2  * 


Vk  "  Vkl  +  Jvk2 


(a) 


(b) 


(12) 


and  define  the  column  vector 


=  (ukl*  Uk2 ’  vkl’  Vk2}' 


(13) 


(the  T  denotes  transpose),  the  output  of  the  kth  quadradic  combiner,  d  , 
can  be  written  as 


dk  =  ^k  Q  ^k 


(14) 


where  Q  has  the  form 


Q  = 


a  0 

0  a 


c ,  -c  b 

1  2 


C2  C1 


C2 

Cl 

0 

b 


(15) 
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The  decision  variable  D  then  is 


D 


=  S  -k  Q 


k-1 


(16) 


in  which  the  L  terms  are  mutually  independent. 

Defining  terms  appropriately,  the  error  probability  takes  the  form 

Pg  =  ProbjD  <  O)s^(t)  transmitted^  .  (17) 

To  make  the  problem  tractable,  we  take  advantage  of  the  slow  variation  of 
the  focus  component  and  assume  that  exp  i(i(t)  is  constant  over  the  sample 
interval  0  £  t  £  T.  With  this  assumption,  the  joint  distribution  of  the 
random  vector  xk  is  Gaussian,  conditioned  on  the  value  of  the  log-normal 
part.  We  then  determine  that 

peU  =  Prob  {D  <  0  I  si(t),  e't'k(t)  =  >  U8) 

and  average  over  a  to  ev-’.luate  the  unconditioned  probability  P  . 

To  determine  the  probability  distribution  of  D,  we  first  compute 

T 

the  characteristic  function  for  d,  =  x,  Q  x  .  This  is  a  quadratic  form 

k  -k  -k 

in  a  Gaussian  random  vector;  the  characteristic  function  is  known  and 
given  by 


cd(0 


exp  j 

I  -  (I  -  2j£RkQ)"1 

*k} 

det 

{i  -  uciyjj 

i 

(19) 


where  is  the  covariance  matrix 


Rk  =  (^k  “  ^k)(4  •  *k> 


(20) 
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Equation  (19)  can  be  written  in  a  form  more  suitable  for  computation 
in  terms  of  the  eigenvalues  and  eigenvectors  of  the  matrix  RQ.  Using 
simple  matrix  manipulations, 

T  -1— 

-x  QMr  M  x 

Cd^)=-fJ= - <21> 

n  (i  -  2ja.r^ 

i=l 

where  the  X^  are  the  eigenvalues  of  RQ,  and  M  is  the  matrix  of  eigen¬ 
vectors  of  RQ  (sometimes  called  the  modal  matrix).  That  is: 

RQ  =  M  A  M-1  ,  (22) 

where 

A  =  diagjx^  ;  (23) 


also, 


T  =  diag< 


-  2JCA.I 


(24) 


The  relationship  between  the  characteristic  function  and  the  error 
probability  is  given  by  Parzen  (1960): 


P 

e 


♦ 


(25) 


where  CQ(£)  is  the  characteristic  function  of  the  decision  variable  D. 
Using  the  fact  that  the  diversity  outputs  are  independent  and  identically 
distributed,  we  have 


P 

e 


I  ♦ 


(26) 
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Averaging  over  f  (and  interchanging  orders  of  integration  and  averaging) 


we  obtain 


(27) 


♦ 
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Ill  NUMERICAL  EVALUATION  OF  P 

e 


A.  Introduction. 

Equation  (27)  is  the  complete  general  expression  for  the  error  rate 
using  the  two-component  channel  model.  Like  many  general  expressions, 
its  evaluation  is  difficult  and  costly.  The  reason  in  this  case  is  that 
Equation  (27)  requires  performing  a  three-dimensional  numerical  integra¬ 
tion  to  a  very  high  degree  of  accuracy,  as  is  given  by  the  difference 
between  one-half  and  the  result  of  the  integration.  Although  we  have 
performed  this  integration  for  several  examples,  we  also  have  considered 
several,  more  easily  evaluated,  special  cases  that  provide  insight  into 
the  relationship  between  P^  and  the  channel  scintillation  parameters. 
These  special  cases  are  obtained  by  turning  off  one  of  the  two  portions 
of  the  model.  The  resulting  two  single-component-channel  models--the 
scatter  (only)  channel  shown  in  Figure  7  and  the  focus  (only  channel 
shown  in  Figure  8--are  important  in  their  own  right.  The  scatter  channel 
can  characterize  Rayleigh  or  Rician  fading  channels  as  well  as  more 

general  situations  where  o  £  a  and  p  /=  0.  The  focus  (only)  channel 

xy  xy 

is  the  model  for  log-normal  fading. 


h  (t) 
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FIGURE  7  SCATTER  (only)  CHANNEL  MODEL 


FIGURE  8  FOCUS  (only)  CHANNEL  MODEL 

In  the  following  sections,  we  present  results  of  our  numerical 
evaluations  of  for  both  single-component  channels  and  the  two-component 
channel  for  FSK  and  DPSK  communication  systems.  Some  of  the  details  of 
the  receiver  setup  (i.e.,  the  way  in  which  the  receiver  was  specialized 
for  the  two  types  of  modulations)  and  the  computation  procedures  are 
given  in  Appendix  A.  For  all  computations,  the  noise  was  assumed  to  be 
stationary  white  Gaussian. 

B.  Single-Component  Channels 

1.  the  Scatter  Component 

For  the  scatter  channel,  the  received  signal  is  given  by 


r(t)  =  [T]  +  hg(t)]  s(t)  +  n(t ) 


This  model  is  perhaps  the  most  useful  of  the  two  because  it  reduces  to 

2  2 

the  Rician  channel  when  a  =  a  and  to  the  Rayleigh  channel  with  the 

x  y 

added  constraint  T]  =  0 .  Furthermore,  T]  =  1  corresponds  to  the  non¬ 
fading,  additive  noise  channel. 

To  characterize  the  received  signal  r(t)  we  must  know  the  following 
parameters--in  addition  to  the  input  waveform  s(t): 

2 

•  T|  (the  energy  in  the  coherent  por tion--note  that  this  fixes 
the  energy  in  the  scatter  portion--see  Eq.  4). 

•  The  SNR. 


a 

— j  (the  ratio  of  the  in-phase  power  to  the  quadrature  power). 
a 

y 
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•  p  (the  correlation  coefficient  between  the  in-phase  and 
quadrature  terms „ 

•  t/T  [the  correlation  time  of  the  scintillation  (t)  relative 
to  the  bit  duration  (T) ;  t-1  is  the  half-power  bandwidth 
for  a  gaussian  fading  spectrum  (see  Appendix)]. 


a.  Coherent  Component  Variations 

The  lossless  channel  [Eqs.  (7),  (8)]  has,  in  effect,  normalized  the 

channel  so  the  coherent  component  T)  is  restricted  to  the  range  0  to  1. 

2 

Figure  9  shows  how  the  power  in  the  coherent  component,  Tj  ,  affects 
at  an  SNR  of  20  dB  for  slow  and  fast  fading.  Observe  that,  for  slow 

fading,  the  error  rate  for  DPSK  is  always  lower  than  the  FSK  error  rate. 

2  ~ 

For  fast  fading,  however,  the  FSK  error  rate  is  lower  for  T)  <  0.6.  The 

reason  for  this  behavior  is  that  FSK  is  inherently  better  than  DPSK  in 

fast-fading  regimes  (Bello  and  Nelin,  1962b),  but  DPSK  is  better  when  the 

fading  is  slow  or  non-existent.  As  the  magnitude  of  the  coherent  component 

increases,  the  effects  of  fading  become  less  and  less  significant,  and  DPSK 

becomes  the  lower  error-rate  system.  Note  also  that  P^  does  not  begin  to 

show  a  significant  reduction  until  over  one-half  the  power  is  in  the  co- 
2 

herent  component  (Tj  >  0,5).  This  is  of  significance  for  obtaining  ap¬ 
proximate  error  rates  for  situations  in  which  some  coherent  component  is 
present,  but  the  fading  dominates. 


2  2 

o  /ct  Variations 
x  y 


The  division  of  the  scintillating  portion  of  the  received  signal 

2  2 

into  in-phase  and  quadrature  components  is  given  by  the  ratio  u^/a^. 

If  the  scatter  ellipse  (see  Figure  3)  is  oriented  vertically,  i.e. 

p  =0,  then  this  ratio  is  also  equal  to  the  axial  ratio  of  the  ellipse, 
xy  r 

For  weak  scintillation,  the  in-phase  component  corresponds  (in  an  approxi¬ 
mate  sense)  to  amplitude  scintillation,  and  the  quadrature  component  cor¬ 
responds  to  phase  scintillation.  Figure  10  shows  the  relationship  between 

2  2  2 
P  and  a  /e  for  slow  and  fast  scintillation  for  T)  =  0.75,  p  =0  and 
e  x  y  xy 
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a  20  dB  SNR.  The  error  rate  drops  rapidly  as  signal  power  is  transferred 
from  the  in-phase  to  the  quadrature  component  for  both  forms  of  modulation 
and  both  scintillation  rates.  This  figure  gives  an  indication  of  the  harm¬ 
ful  effect  of  amplitude  fading  relative  to  phase  scintillations--at  least 
for  FSK  and  DPSK.  Coherent  phase  systems  (coherent  PSK.) ,  which  are  known 

to  be  more  sensitive  to  rapid  fading,  would  probably  not  show  such  a 

2  2 

dramatic  reduction  in  P  as  a  /a  decreases.  Figure  11  shows  P  as  a 

e  2  x  y  °  e 

function  of  the  SNR  with  a  /a^  as  a  parameter. 

x  y 


c.  Correlation  Coefficient  Variations 

The  correlation  coefficient  p  has  significant  effect  on  the  error 

xy 

rate  as  shown  in  Figures  12  and  13.  (The  severe  drop  in  P  occurs  for 

e 

larger  values  of  p  (>  0.5)  and  is  probably  not  a  situation  likely  to 

occur  in  real  life.)  The  rapid  fall  of  P  as  p  ->  1  is  because  the 

e  xy 

random  scintillation  is  becoming  more  ordered  and,  in  the  process,  de¬ 
veloping  a  lower  limit  on  the  possible  amplitude  excursions.  For  example, 

the  line  L  in  Figure  14  is  the  limit  of  the  scatter  plot  as  p  —  -1  when 
2  2  2 

7]  =0.75  and  o^/a^  =0.5.  The  signal  amplitude  can  never  drop  below  the 

distance  of  closest  approach  of  L  to  the  origin  (the  line  OA) ;  hence  the 
amplitude  fading  has  developed  a  lower  limit.  This  lower  limit  is  im¬ 
portant  because  the  main  reason  fading  increases  the  error  rates  as  much 
as  it  does  is  the  finite  probability  that  a  deep  fade  will  occur  and  the 
signal  will  drop  into  the  noise. 


d .  Fading  Rate 

The  influence  of  fading  rate  on  P^  has  been  examined  for  Rayleigh 

fading  by  Bello  and  Nelin  (1962b);  however,  we  have  been  unable  to  find 

a  study  of  the  fading  rate  in  connection  with  any  other  type  of  fading 

(such  as  Rician).  Bello  and  Nelin  show  that  the  fading  rate  establishes 

a  lower  limit  to  the  error  rate,  termed  the  irreducible  probability  of 

error  (P  ).  Once  this  limit  is  reached,  further  increases  in  the  SNR 
IK 

will  not  result  in  reductions  in  P  . 

e 
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FIGURE  14  SCHEMATIC  DRAWING  OF  SCATTER  DIAGRAM 
ELLIPSE  AS  Pv„  -*  -1 

An  irreducible  probability  of  error  is  also  established  by  the 

2 

scintillation  rate  of  our  scatter'  channel,  given  fixed  values  of  T]  , 

2  2 

p  and  ct  /o  .  Figure  15  shows  the  PTD  versus  the  relative  fading 
xy  x  y  IR  2  2 

bandwidth  (T/t)  for  Rician  conditions  (a  /a  =  1,  p  =0)  for  several 
2  2  x  y  xy 

values  of  T]  including  T]  =0  (Rayleigh  fading).  It  is  interesting  that 

the  Rician  curves  show  the  same  slope  as  the  Rayleigh  curves,  differing 

only  by  a  fixed  offset  due  to  the  coherent  component. 

Figure  16  shows  P  for  non-Rician  conditions.  The  behavior  is 
IK 

much  the  same  as  the  Rician  channel  shown  in  Figure  15  except  for 

2  2 

ct  /a  =  0.2.  Here,  the  difference  between  DPSK  and  FSK  increases 
x  y 

dramatically  as  the  fading  bandwidth  increases  with  the  DPSK  system 
obviously  affected  by  the  rapid  phase  changes  due  to  the  scintillation. 
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The  flattening-out  at  the  low  end  of  the  curves  is  probably  because  this 
range  of  is  beyond  the  limit  of  accuracy  of  the  integration. 


The  Focus  Component 


The  other  half  of  the  two-component  model,  the  focus  channel,  is 

essentially  a  slow-fading  log-normal  channel.  We  have  only  performed 

computations  for  slow  log-normal  amplitude  fading  (i.e.,  a  £  0, 

X 

a  =  0,  p  =0),  shown  in  Figure  17.  These  are  identical  with  results 
9  X<P 

obtained  by  numerical  integration  of  the  standard  expression  for  the 
probability  of  error  in  a  slow-fading  channel 


P 

e 


f(x)  P(x)dx 


(29) 


where  x  is  the  SNR,  P(x)  is  the  log-normal  probability  distribution, 
f(x)  is  the  probability  of  error  for  DPSK  or  FSK  in  a  non-fading  channel, 
i.e., 


Two- Component  Channels 


The  error  rate  for  the  two-component  channel  was  determined  for  the 
examples  of  moderate  and  severe  scintillation  used  by  Scott  and  Knepp 
(1978)  of  MRC.  The  channel  parameters  and  resulting  error  rates  are 
given  in  Table  1.  The  error  rates  we  obtain  for  moderate  scintillation 
are  much  higher  than  those  obtained  by  MRC  and  cannot  be  explained  by 
the  difference  in  modulation  schemes.  Indeed,  DPSK  or  FSK  should  perform 
better  than  PSK  in  the  rapid  scintillation  of  this  channel.  A  possible 
cause  of  this  discrepancy  is  the  small  sample  size  used  by  Scott  and  Knepp 
(they  used  a  Monte  Carlo  simulation  with  a  sample  size  of  ~  3000  to  esti¬ 
mate  P  .  Our  approach  does  not  have  this  limitation, 
e 
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Table  1 


SUMMARY  OF  CHANNEL  PARAMETERS 


Channel  Parameter 

Moderate  Scintillation 

Severe  Scintillation 

MRC* 

SRI* 

MRC* 

SRI1- 

Scatter  Component 

S4  index 

0.5920 

0.9825 

Tl2 

0.411 

0.0174 

CTx 

0.2922 

0.9164 

CTy 

0.4585 

0.8756 

CTx/ay 

0.4061 

0.4061 

1.0952 

1.0952 

pxy 

-0.2039 

-0.2039 

-0.0022 

-0.0022 

T/T* 

0.041 

0.  Ill 

Focus  Component 

S^  index 

0.2569 

0.4416 

ax 

0. 1282 

0. 1282 

0.209 

0.209 

a® 

1.3328 

1.3328 

3.6563 

3.6563 

px* 

-0.5435 

-0.5435 

-0.1376 

-0.1376 

*/T* 

0.0013 

0 

0 

Error  Probabilities 

P  (PSK)$ 
e 

0.39  x  10‘3 

0.23  x  10"1 

-1 

Pg  (DPSK) s 

0.197  X  10 

0.865  X  10 

P  (FSK)5 

0.151  x  10'1 

0.345  X  10- 

e 

* 

MRC  =  Values  given  by  Scott  and  Knepp  (1978). 

SRI  =  Normalized  MRC  values  suitable  for  our  computations. 
*Scott  and  Knepp  used  a  data  rate  of  75  bits/sec. 

R 

SA  mean  SNR  of  16.25  dB  was  used. 


PROBABILITY  OF  ERROR 


In  addition  to  the  computations  shown  in  Tuble  I,  we  varied  the 

focus  parameters  (a  ,  a  and  p  )  above  and  below  the  moderate  scintilla 

X  $ 

tion  levels  and  computed  the  resulting  error  rates.  Neither  a  nor  p 

had  any  effect  whatsoever  on  P  .  This  is  understandable  when  one  recall 

e 

that  DPSK  is  a  differential  system  and  should  be  transparent  to  slow 
phase  variations.  Of  course,  FSK  will  also  be  transparent  to  slow  phase 
variations.  On  the  other  hand,  a  system  employing  fixed-reference  PSK 
will  be  severely  degraded  by  any  phase  variations--f ast  or  slow. 

The  error  rate  was  affected  by  variations  of  a  ,  as  shown  in 

X 

Figure  18.  However,  P^  is  not  a  sensitive  function  of  this  parameter 
at  the  level  of  the  other  scintillating  parameters  of  this  example. 

(The  principal  determinants  of  P^  in  this  example  are  the  SNR  and  the 
fading  rate.) 


IV  CONCLUSIONS 


The  two- component  channel  model  provides  both  a  means  of  calculating 
the  probability  of  error  for  the  transionospheric  channel  and  an  insight 
to  the  mechanism  of  error  generation.  Although  extensive  evaluation  of 
the  complete  expression  for  P  (to  study  a  system  under  many  different 
conditions,  for  example)  would  be  costly,  it  is  feasible.  Also,  the 
model  lends  itself  to  simplifications  that  allow  special  cases  to  be 
studied  economically  and  easily. 

Real  channels  do  not  have  an  unlimited  range  of  possible  parameter 
variations.  Furthermore,  the  statistics  of  these  parameters  (means, 
variances,  etc.)  are  usually  not  independent  of  one  another  even  if  the 

parameter  values  themselves  are  independent.  For  example,  a  high  scatter 

2  2 
scintillation  index  (low  7]  )  will  usually  mean  a  high  value  of  a  as 

well.  Hence,  studying  a  system  under  conditions  that  are  likely  to  occur 

in  operation  does  not  necessarily  mean  performing  an  inordinately  large 

number  of  computations. 

The  dominant  source  of  errors  induced  by  transmitting  through  the 
two-component  channel  is  the  scatter  component.  The  scatter  component 
controls  the  amount  of  coherent  signal  and  the  amplitude  and  rapidity  of 
the  fast  scintillations.  These  parameters,  together  with  the  noise,  are 
the  most  critical  aspects  of  the  signal.  The  focus  component  has  only  a 
secondary  effect  on  P  ,  at  least  when  it  is  present  in  the  same  relative 
proportions  as  seen  in  our  examples.  For  the  modulation  schemes  we  have 
examined,  the  focus  term  can  almost  be  viewed  as  a  simple  log-normal 
perturbation  on  the  scatter  component. 

Although  Eq.  (27)  is  expressed  in  terms  of  an  L-branch  diversity 
receiver,  all  of  our  analysis  has  been  for  L  =  1,  i.e.,  no  diversity. 

This  is  because  considerable  study  has  already  been  directed  toward  the 
effects  of  diversity  reception  and  the  general  behavior  of  diversity 
systems  is  reasonably  well  understood.  Our  principal  intention  has  been 
to  study  the  way  in  which  various  channel  parameters  affect  the  error  rate. 
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Appendix 


DETAILS  OF  THE  ERROR-RATE  CALCULATION 


FSK  Receiver 

For  FSK,  Che  information  and  reference  filters  in  Figure  2  are 
matched  to  the  mark  and  space  waveforms;  hence,  the  impulse  responses  are 


-jrntt  v 

gl(t)  =  sL  (t)  =  V2E/T'  e  T  ) 

>  0<t<T  (A-l) 

*  _  +jnnt  \ 

gQ(t)  =  sQ  (t)  =  JlETf  e  T  / 


where  E  is  the  total  energy  in  the  pulse,  T  is  the  pulse  duration,  and 
n  is  an  integer  equal  to  the  frequency  separation  between  the  mark  and 
space  frequencies  normalized  with  respect  to  1/T  (the  asterisk  denotes 
the  complex  conjugate).  The  quadratic  combiner  matrix  is 


Q  = 


10  0  0 

0  10  0 

0  0-10 
0  0  0  -1 


(A-2) 


DPSK  Receiver 

For  DPSK  the  impulse  responses  of  the  receiver  filters  are 

gl(t)  =  s  "(t)  -  s0"(t)  =  V2E/T'  0<t<T  (A-3) 


gQ(t)  =  gL(t  +  T) 


-T<t<0  ,  (A-4) 
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and 


s,(t)  -„'2E/f 


0<t<T 


sQ(t)  =  -v'2E/T 


(A- 5) 


The  quadradic  combiner  matrix  is 


Q  = 


0  0  10 

0  0  0  1 

10  0  0 
0  L  0  0 


(A- 6) 


Hence, 


and 


"  fh (t 


)r  (t)dt 
k 


o 

■  fh 


(t  +  T)rk(t)dt  = 


A 


(t)rk(t  -  T)dt 


(A- 7) 


(A-8) 


Note  that  the  reference  filter  is  identical  to  the  information  filter 
with  the  input  delayed  one  bit  duration. 


R  Matrix  Computation 

It  is  necessary  to  compute  the  covariance  matrix 


R  =  (x  -  x)  (x  -  x)"^  .  (A- 9) 

This  matrix  is  composed  of  auto-  and  cross-correlations  of  the  real  and 
imaginary  parts  of  u  and  v  (see  Eqs.  9  and  10)  such  as 
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r  11  =  (U1  -  ulr 


ri2  =  r21  =  (U1  ~  G1)(U2  -  u2> 


(A-10) 


r 13  =  r 31  =  (U1  '  U1)(V1  -  V 


et  cetera. 


The  R  matrix  is  completely  determined  once  the  following  six  relations 
are  known: 


(u  - 

u)2  = 

,  2  -2 
(U1  “  Ul} 

-  <U2 

-  u2)  +  2j (u^u2  - 

•  V2} 

(A-lla) 

, 

, 

2  -2 

~2 

-2, 

(u  - 

u)(u"'  - 

u  )  = 

(U1  -  V 

+  (u2 

-  u2) 

(A- lib) 

oT 

-  u)  (v 

~  = 

(uTi  -  r‘ 

l"l>  ' 

~2 

-  u,vj  + 

j(V2 

-  U  V  ) 

i  r 

+  j(u2Vl 

-  vV 

(A-llc) 

(u  - 

U)(v*  - 

v“)  = 

(ulVl  -  u 

iV  + 

(U2V2 

'  VV  - 

J(U1V2 

-  VP 

+  j(u2Vl 

-  vV 

(A-lld) 

.  2 

,~2  -2 

~2 

-2, 

(V  - 

v)  = 

(vi  -  V 

-  (V2 

-  V 

+  j(Vlv2  - 

■  VP 

(A-lle) 

(V  - 

v)(v"  - 

V*)  = 

~2  —2 
(v2  -  vx) 

+  (V2 

V 

, 

(A-llf) 

Note  that  any  element  of  R  can  be  obtained  by  taking  the  real  or  imaginary 
part  of  a  sum  or  difference  of  two  of  the  above  expressions.  For 
example, 
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r23  =  (V1  •  V1)(U2  '  V  =  V1V2  -  V1U2 


=  1/2  Im|(u  -  u)(v  -  v)  +  (u  -  u)(v  -  v  )| 


(A- 12) 


For  the  averages,  we  have 


where 


u  =y'gl(t)r(t)dt  =  y'g1(t)r(t)dt  , 


r(t)  =  [T|  +  h(t) ]  e^t^s(t)  +  n(t) 


(A-13) 


=  7]  e^(t)s(t) 


(A- 14) 


This  average  is  intended  to  be  taken  over  a  period  of  time  that  is  long 
with  respect  to  a  bit  duration  T,  but  short  with  respect  to  the  time  scale 
of  the  focus  component.  Hence,  is  a  constant  over  the  averaging 

period  and 


r (t )  =  Tle^0  s(t) 


(A-15) 


A  similar  result  is  obtained  for  v. 


Equations  (A-ll)  are  easily  evaluated  in  terms  of  the  channel  corre¬ 
lation  functions  h(t^)h(t2)  and  h( t^ )h*(tp . 
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For  example 


(u  -  ur  = 


■L 

/ dtg^CtJ^rCt)  -  T]  eV^t')s(t); 


L  o 

T  T 


/-./  dt2  gl(tl)gl(t2^r^tl^r^t2^ 


o  o 


Lf- 


(t)s(t) 


(A-16) 


Because  \Jf ( t )  varies  slowly,  we  have 


r(t1)r(t2)  =  |  [T]  +  h(t^)  ]ev^t^s(t^>  +  n  ( t  x  )|  j[T]  +  h(  t2)  ]e^ t  ^  s  ( t2>  +  n(t2>^ 

=  [ll2  +  hU1)h(t2)  J  e  ^ v  ( 1 1 }  +  ’;,(t2)]s(t1)s(t2)  +  nU1)a(t2) 


=  e 


2*(t)  [l]2  +  h(tl)h(t2)  s(t1)s(t2)  +  n(t1)n(c2)  .  (A- 17 ) 


The  noise  is  a  complex  process--the  real  and  imaginary  parts  of  which  are 
(1)  independent,  (2)  zero  mean,  and  (3)  identically  distributed;  therefore 


n(tl)n(t2)  =  [x^t^  +  jyn(t1)^xn(t2->  +  jyn(t2)' 


=  Xn(tl)xn(t2)  '  yn(tl)yn(t2)  =  ° 


(A- 18) 


Hence, 


-  2  2 
(u  -  u)  =  e 


fdtl  f  dt2  gl(tl^gi<'t2^Ahh^tl’t:2^S^tl^S^t:2^  ’ 

Jo  Jo  (A- 19) 
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where  A^h( c ^ ^  =  h(t^)h(t2).  The  other  expressions  in  (A-ll)  are  de¬ 
termined  similarly  (the  noise  term  obviously  does  not  always  cancel  it¬ 
self  out). 

The  channel  correlation  functions  are  complex  and  given  by 


h(t1)h(t2)  =  [x(tL)  +  jy(t1)][x(t2  +  jy(t2)] 

=  x(t1)x(t2)  -  y(tx)y(t2)  +jjx(t1)y(t2)  +  xC^yCt^j  ^A_20a) 


h( c  j_ )  h*(  1 7)  =  x(t1)x(t2)  +  y(t1)y(t2)  .  (A-20b) 

I' he  channel  is  assumed  to  be  stationary,  hence  Eqs.  (20a, b)  are  functions 
ot  only  the  difference  T  =  t ^  -  t2>  e.g., 


x(t  )x(t2)  =  r3QC(^t> 


(A-21) 


We  also  assume  that  all  channel  correlation  functions  are  Gaussian  of  the 
form 


2  2 

R(*t)  =  <r2e 


(A-22) 


where  the  constant  a  is  given  by 


2  2 
o  =  x^(t) 

XX 


2  2 
o  =  y  (t) 

yy  7 


a  =  p  a  a 
xy  xx  yy 


and 


p  =  x(t)y(t) 
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This  particular  form  is  taken  from  Belo  and  Nelin  (1962).  We  note 
that  B  =  1/t  is  the  half- power  bandwidth  of  the  fading  channel;  thus> 
is  not  the  usual  e-^  point  of  the  gaussian  correlation  function. 
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